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Confirmation of solar-like oscillations in eta Bootis 

H. Kjeldsen, 1 ' 2 T. R. Bedding, 3 I. K. Baldry, 3 ' 4 H. Bruntt, 2 R. P. Butler, 5 D. A. Fischer, 6 
S. Frandsen, 2 E. L. Gates, 6 F. Grundahl, 2 K. Lang, 1 G. W. Marcy, 6 A. Misch, 7 and S. S. Vogt 7 

ABSTRACT 

We obtained time-series spectroscopy of the GO subgiant r] Boo in an attempt to 
confirm the solar-like oscillations reported by Kjeldsen et al. (1995). We recorded 1843 
spectra over six consecutive nights with the Nordic Optical Telescope, which we used 
to measure equivalent widths of strong temperature-sensitive lines. We also measured 
velocities from 1989 spectra obtained through an iodine reference cell at Lick Obser- 
vatory over 56 nights that were badly affected by weather. Our analysis also included 
velocity measurements published by Brown et al. (1997) and the original equivalent- 
width measurements by Kjeldsen et al. (1995). All four data sets show power excesses 
consistent with oscillations, although with a range of amplitudes that presumably re- 
flects the stochastic nature of the excitation. The highest peaks show regularity with 
a large separation of Av = 40.4 /^Hz and we identify 21 oscillation frequencies from the 
combined data. 

Subject headings: stars: individual {rj Boo) — stars: oscillations — techniques: radial 
velocities 

1. Introduction 

The search for solar-like oscillations is finally yielding success. Observations of the sub-giants 
Procyon (Martic et al. 1999; Barban et al. 1999) and (5 Hydri (Bedding et al. 2001; Carrier et al. 
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2001) have shown very good evidence for oscillations. More recently, there has been an unam- 
biguous detection of p-mode oscillations in the main sequence star a Cen A by Bouchy & Carrier 
(2001, 2002). All these results were based on velocity measurements obtained using high-dispersion 
spectrographs with stable reference sources. 

Another method for detecting oscillations was suggested by Kjeldsen et al. (1995, hereafter 
Paper I) . This involved monitoring changes in the equivalent widths (EWs) of temperature-sensitive 
spectral lines. In Paper I we reported evidence for oscillations in the GO sub-giant r\ Boo, based 
on measurements of Balmer-line EWs. We presented this as the first clear evidence of solar-like 
oscillations in a star other than the Sun. The observations were obtained over six nights with the 
2.5 m Nordic Optical Telescope on La Palma, and consisted of 12684 low-dispersion spectra. In the 
power spectrum of the equivalent-width measurements, we found an excess of power at frequencies 
around 850 fiHz. The average amplitude inferred for the oscillations was about 7 times greater 
than solar, in rough agreement with the empirical scaling relation suggested by Kjeldsen & Bedding 
(1995). Comb analysis of the power spectrum, described in Paper I, suggested a regular spacing of 
Au = 40.3 //Hz. Based on this, we identified thirteen oscillation modes. Similar observations of the 
daytime sky showed the five-minute solar oscillations at the expected frequencies. 

The frequencies for rj Boo reported in Paper I, taken with available estimates of the stellar 
parameters, were in good agreement with theoretical models (Christensen-Dalsgaard et al. 1995; 
Guenther & Demarque 1996). Particularly exciting was the occurrence in theoretical models - and 
apparently in the observations - of 'avoided crossings,' in which mode frequencies are shifted from 
their usual regular spacing by effects of gravity modes in the stellar core. Since then, the improved 
luminosity estimate for 7] Boo from Hipparcos measurements has given even better agreement with 
the measured value of Av (Bedding et al. 1998). 

Meanwhile, a search for velocity oscillations in rj Boo by Brown et al. (1997) has failed to 
detect a signal, setting limits at a level below the value expected on the basis of the EW results. 
Although the data were sparse (22 hours spread over 7 successive nights) and the precision was 
degraded by the relatively fast rotation of the star (v sin/ = 13kms _1 ), the analysis by Brown et al. 
(1997) was careful and thorough, and the results seem to be inconsistent with those of Paper I. 
More recently, Carrier et al. (2003, and paper in preparation) reported velocity measurements using 
the CORALIE and ELODIE spectrographs that showed a clear excess of power and a frequency 
spacing of 39.6 //Hz. 

In this paper we present additional observations of rj Boo, obtained in 1998 in both EW and 
velocity. We also re-analyse our 1994 EW measurements and the velocity measurements of Brown 
et al. (1997). We confirm our earlier claim for oscillations in r\ Boo and identify more than twenty 
oscillation frequencies from the combined data. 
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2. Data 

2.1. Equivalent- width observations (NOT98) 

We observed rj Boo over six nights during May 1998 using ALFOSC (Andalucia Faint Object 
Spectrograph and Camera) on the 2.5 m Nordic Optical Telescope on La Palma. This is the same 
telescope used in Paper I but with a different spectrograph. We obtained seven echelle orders 
covering the range 370-700 nm at a dispersion of 0.04 nm/pixel. The CCD was a Loral 2k by 2k 
device, of which we used 700 by 1200 pixels. 

Spectra were taken with a typical exposure time of lis and a dead time between exposures 
of 17 s. They were averaged in groups of three before writing to disk, resulting in a total of 1843 
spectra (sampling rate 1/84 s) in 44.2 hours over six consecutive nights (1998 May 1-6). The 
distribution of spectra over the six nights was: 88, 222, 377, 391, 358 and 407. 

We used the method described by Kjeldsen et al. (1999) to measure equivalent widths of six 
strong temperature-sensitive lines: Ha, H/3, FLy, Mgi, Nai and Fei. A weighted mean of these six 
values was calculated, taking into account the differing temperature sensitivities of the lines. The 
resulting time series is shown in Fig. 1. 

2.2. Velocity observations (Lick98) 

We also observed rj Boo in 1998 with the Hamilton Echelle Spectrometer and the 0.6-m Coude 
Auxiliary Telescope (CAT) at Lick Observatory (Vogt 1987). To produce high-precision velocity 
measurements, the star was observed through an iodine absorption cell mounted directly in the 
telescope beam. 

We were allocated 56 of the 59 nights from 1998 April 6 to June 3, but the weather was 
unseasonably poor, permitting observations on only 26 nights (and many of these were partly lost). 
The exposure time was 120 s, with a dead time between exposures of the same amount. On the 11 
best nights we obtained 95-120 spectra per night (sampling rate: 1/245 s), and the total number 
of spectra obtained was 1989 (about one third of that possible with no weather losses). 

Extraction of radial velocities from the echelle spectra followed the method described by Butler 
et al. (1996). As mentioned in the Introduction, the precision is degraded by the relatively fast 
stellar rotation. The star is a spectroscopic binary with a period of 494 d (Bertiau 1957), and the 
orbital motion was removed from the velocity time series by fitting and subtracting a fifth-order 
polynomial. The resulting velocity measurements are shown in Fig. 2. We are confident that 
the long-term velocity variations are not instrumental, given that velocities for r Ceti, which we 
observed on most of the nights, were stable at the 5ms _1 level. These night-to-night variations in 
rj Boo are presumably due to stellar activity, which is commonly observed in rotating G-type stars 
(see, e.g., Saar et al. 1998; Santos et al. 2000). 
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2.3. Published equivalent-width observations (NOT94) 

We have included in this analysis the time series of 12684 EW measurements obtained with 
the Nordic Optical Telescope in 1994. These are identical to the data presented in Paper I, with 
the exception that a high-pass filter was not applied. The result is an increase in noise at low 
frequencies, as expected from a 1// noise source, which more accurately reflects the actual stellar 
and instrumental noise (for more details, see Bedding & Kjeldsen 1995). 

2.4. Published velocity observations (AFOE95) 

We have also analysed 555 velocity measurements of r\ Boo obtained with the AFOE spectro- 
graph during 22 hours spread over seven successive nights in 1995 March. These measurements 
were described by Brown et al. (1997) and were kindly provided to us in electronic form by Tim 
Brown. 

3. Analysis and discussion 

The power spectrum of each time series was calculated as a weighted least-squares fit of si- 
nusoids (Frandsen et al. 1995; Arentoft et al. 1998), with a weight being assigned to each point 
according its uncertainty estimate. The results for the four data sets are shown in Figs. 3-6. In 
each case, we show both the conventional power spectrum (upper panels) and a smoothed version 
in which the vertical scale has been converted to power density (lower panels). As discussed by 
Kjeldsen & Bedding (1995, Appendix A.l) and Kjeldsen et al. (1999, Sec. 5), the conversion to 
power density is achieved by multiplying by the effective observing time, which we calculated in 
each case by integrating under the spectral window. 

Some level of excess power in the range 600-1100 /iHz is apparent in all four data sets, although 
it is strongest in the NOT94 data. Further discussion of the measured power levels is given below 
in Sec. 3.4. First, however, we will discuss the oscillation frequencies. 

Mode frequencies for low-degree solar-like oscillations are approximated reasonably well by the 
asymptotic relation: 

!/„,, = Av(n+\l + e)-l(l + l)D . (1) 

Here, n and I are integers that define the radial order and angular degree of the mode, respectively; 
Au (the so-called large separation) reflects the average stellar density, D is sensitive to the sound 
speed near the core and e is sensitive to the surface layers. It is conventional to define 5vq2, the so- 
called small separation, as the frequency spacing between adjacent modes with / = and / = 2. We 
can further define Svqi to be the amount by which I = 1 modes are offset from the midpoint between 
the I = modes on either side. If the asymptotic relation holds exactly, then it is straightforward 
to show that Dq = ^Si>q2 = 5^01 • 
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3.1. Extraction of frequencies 

We extracted the frequencies of the strongest peaks in each power spectrum in the range 600- 
1100 juHz. We used a simple iterative algorithm, in which the highest peak was identified and the 
corresponding sinusoidal variation was subtracted from the time series. The power spectrum of the 
residuals was then calculated and the process was iterated until all peaks with amplitudes more 
than 2.5 times the noise floor had been extracted. The numbers of peaks extracted from the four 
data sets are summarized in Table 1. 

Setting the threshold at 2.5a gives us a chance to detect the weaker oscillation modes, but it 
also means we will select some noise peaks. To investigate this, we have conducted simulations in 
which we analysed noise spectra that contained no signal. The last column in Table 1 shows the 
number of peaks found above the 2.5a threshold. This indicates that about 6 of the 35 detected 
peaks are expected to be due to noise. 

3.2. Large frequency separation 

We next investigated whether the extracted frequencies have a regular spacing, as is expected 
for p-mode oscillations. In Paper I, we described a comb analysis of the NOT94 power spectrum 
that revealed a regular spacing of Av = 40.3 /iHz. Here, we analyse the 22 frequencies from the 
other three data sets (AFOE95, NOT98, Lick98) using autocorrelation, as follows. Firstly, each 
extracted frequency was allocated a power corresponding to the square of its S/N, in order to 
give higher weight to the more reliable peaks. Secondly, to allow for the likelihood that some 
of the extracted frequencies should be shifted by ±ld _1 (±11.57 fiHz), we extended the table of 
frequencies by a factor of three by including these side-lobes (but with half the power of the central 
peaks). We then calculated the autocorrelation of these 66 frequencies. This is shown, smoothed 
to a resolution of 3.5 /xHz, as the solid line in Fig. 7. 

The three peaks in the autocorrelation correspond to the large separation and its daily aliases, 
leading us to estimate a value of Av = 40.4 /iHz. For comparison, the dotted line in Fig. 7 
shows the autocorrelation for the 13 frequencies extracted from the NOT94 data, which yields a 
large separation of Au = 40.5 fiHz. The excellent agreement between these independent data sets 
confirms the result of Paper I and also agrees well with the value of 39.6 fiHz reported by Carrier 
et al. (2003). In summary, there can be no doubt that the large frequency separation of r\ Boo 
is about 40 ^Hz, which is in excellent agreement with theoretical models (Christensen-Dalsgaard 
et al. 1995; Guenther & Demarque 1996). 
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3.3. Identification of frequencies 

Given the large separation, we next attempted to identify the individual modes, remembering 
that some of them will need to be shifted by ±11.57 /xHz. Figure 8 shows the 35 frequencies from 
all four data sets, displayed in an echelle diagram. We expect modes with a given I value to 
form vertical ridges, for I = 0, 1 and 2. We were able to achieve this by shifting 15 frequencies 
by ±11.57/iHz, as shown in the figure. The final frequencies are given in Table 2. Note that 
the number of peaks that are classified as noise is consistent with the estimates made above (see 
Table 1). 

A possible problem is that, modulo the large separation, peaks corresponding to [y n i — 
11.57/xHz) for I = 1 are only separated by 1.1 /iHz from peaks corresponding to {y n j, + 11.57/xHz) 
for 1 = 2. It is possible that some of those peaks have been shifted the wrong way and are therefore 
wrongly identified. The relevant peaks are 699.6 fiHz (from Lick98), which we identified as an alias 
of an I = 1 mode, and 822.3 /uHz (from NOT94), which we identified as an alias of an / = 2 mode. 
A similar comment, although to a lesser extent, applies to peaks corresponding to {y n ,l ~ 11-57 /uHz) 
for I = and to {y n ,i + 11-57 juHz) for I = 1, which are only separated by 1.9 /iHz. 

Note that we have identified the peak at 749.3 fiHz as possibly being an I = 1 mode that is 
displaced by an avoided crossing. This peak has the second-highest S/N of all the peaks and is 
therefore not likely to be due to noise or to be an alias. In any case, shifting this peak by ±11.57 / uHz 
would not bring it into agreement with I = or 2. 

Some of the frequencies are detected more than once, as indicated by the ditto marks (" ) in 
Table 2. In those cases, we have combined the measurements into a weighted average. The final 
list of 21 frequencies is given in Table 3, and these are shown as an echelle diagram in Fig. 9. The 
uncertainties in the frequencies reflect the S /N of the relevant peak (or peaks) . 

We next estimated the large separation separately for each value of I, and the results are given 
in the last line of Table 3 (the two modes marked with an asterisk, which may be affected by 
avoided crossings, were excluded from the fit). The weighted average of these three Avi yields the 
value of Av shown in Table 4. The other parameters in that table were calculated by fitting the 
same 19 frequencies to the asymptotic relation. All are consistent with the values given in Paper I 
but have smaller uncertainties, thanks to the larger number of detected frequencies. 

In Fig. 9 we also show the 13 frequencies reported in Paper I. We can see that six were recovered 
in the new analysis, while five were not recovered but lie close to one of the ridges (within typical 
uncertainties), so perhaps can still be taken as reliable detections. Finally, two frequencies (786.2 
and 950.3 /iHz) fall well away from the I = 1 ridge and were not confirmed by the revised analysis. 
It is possible that they represent mixed modes that are shifted by avoided crossings, but further 
observations are needed to confirm this. 
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3.4. Oscillation amplitudes 

The above analysis gives strong evidence that the highest peaks in the four power spectra are 
due to solar-like oscillations. However, we have also noted that the amplitude of the power excess is 
somewhat stronger in the NOT94 observations than in the other three data sets. It seems plausible 
that these differences reflect the stochastic nature of the excitation mechanism, but a more definite 
statement is hampered by our limited knowledge of the amplitudes of oscillations in subgiant stars. 
It is also possible that some contribution to the variations is due to an unidentified systematic error 
in the calibration that was used to convert equivalent- width amplitudes to velocities (see Paper I). 

Taken at face value, the observations indicate that peak oscillation amplitudes in r] Boo are 
typically 3-5 times solar. This conclusion is consistent with the upper limits reported by Brown 
et al. (1997) from the AFOE95 data. 

4. Conclusions 

We have presented new observations of r\ Boo in velocity (Lick98) and equivalent width (NOT98) 
that show some evidence for excess power in the range 600-1100 ^Hz. The velocity measurements 
published by Brown et al. (1997, AFOE95) also show a slight power excess when smoothed. None 
of these signals is as strong as the original equivalent-width observations (Paper I, NOT94), which 
may reflect the stochastic nature of the excitation or may indicate a problem with the calibration 
of the NOT94 equivalent-width estimates. 

We extracted the highest peaks in each data set and used the autocorrelation to search for 
regularity. The three newer data sets (AFOE95, NOT98 and Lick98) combined to give a clear 
autocorrelation signal at a spacing of Au = 40.4 /xHz, giving independent confirmation of the result 
of Paper I. This, combined with the recent work by Carrier et al. (2003), leaves little doubt that 
the large frequency separation of i] Boo is about 40 fiHz. This allowed us to identify 21 frequencies 
in rj Boo which have been compared with theoretical models by Di Mauro et al. (2003). The results 
confirm the claim made in Paper I for the first clear evidence of solar-like oscillations in a star other 
than the Sun. Future observations of r\ Boo, particularly with the MOST spacecraft (Matthews 
et al. 2000), should measure more oscillation modes with greater frequency precision. 

We would be pleased to make the data available on request. Please contact Tim Bedding 
(bedding@physics . usyd . edu . au) . 

We thank Tim Brown for providing the AFOE velocity measurements of rj Boo. This work was 
supported financially by the Australian Research Council (TRB and 1KB), National Science Foun- 
dation Grant AST-9988087 (RPB), SUN Microsystems, and the Danish Natural Science Research 
Council and the Danish National Research Foundation through its establishment of the Theoretical 
Astrophysics Center (HK). 
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Fig. 1. — EW measurements of rj Boo obtained at the NOT 98 (lower panel) and the corresponding 
uncertainties (upper panel). 
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Fig. 2. — Velocity measurements of 7] Boo obtained at Lick (lower panel) and the corresponding 
uncertainties (upper panel). 
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Fig. 3. — Power spectrum of the NOT98 EW measurements of rj Boo. 
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Power spectrum of the Lick98 velocity measurements of rj Boo. 
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Fig. 5. — Power spectrum of the NOT94 EW measurements of rj Boo that were published in Paper I. 
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Fig. 6. — Power spectrum of the AFOE95 measurements of 77 Boo that were published by Brown 
et al. (1997). 




Fig. 7. — Solid line: the autocorrelation of the 22 frequencies, together with their daily aliases, that 
were extracted from AFOE95, NOT98 and Lick98 data sets. Dotted line: the same, but for the 13 
frequencies from the NOT94 data. The peaks at 40 /iHz represent the large frequency separation 
of 77 Boo. 
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Fig. 9. — Echelle diagram showing the 21 frequencies from Table 3, with 2a error bars. For 
comparison, we also show the 13 frequencies reported in Paper I. 



Table 1: Results of the frequency analysis 



Dataset Noise level Extracted peaks Noise peaks 

NOT94 15.1 ppm 13 1.3 ± 0.4 

AFOE95 0.67 ms- 1 5 1.0 ± 0.3 

NOT98 16.3 ppm 3 0.5 ± 0.2 

Lick98 0.41 ms- 1 14 3.3 ± 1.0 



17- 



Table 2: Identification of extracted frequencies 
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4.4 




7£Q Q 


1 = 1 


Lick98 


2.6 


( 00. U — 11.0 ( 


7£Q /i 
= fOo.4 


11 


AFOE95 


2.9 




793 1 


I = 1 


NOT98 


2.8 




806.7 


noise 


AFOE95 


3.0 


798.5 + 11.57 


= 810.1 


I = 2 


Lick98 


2.6 


822.3- 11.57 


= 810.7 


11 


NOT94 


4.3 


801.1 + 11.57 


= 812.7 


I = 


NOT94 


3.1 


825.1 - 11.57 


= 813.5 




NOT94 


3.1 




815.9 


noise 


Lick98 


3.0 




849.9 


I = 2 


NOT94 


2.8 




853.6 


I = 


NOT94 


5.6 


905.8- 11.57 


= 894.2 


I = 


NOT94 


3.2 




955.6 


I = 1 


Lick98 


2.5 


960.1 + 11.57 


= 971.7 


Z = 2 


Lick98 


3.0 


962.0+ 11.57 


= 973.6 


Z = 


NOT94 


3.2 




975.7 




NOT98 


2.9 




1034.3 


I = 1 


Lick98 


2.6 




1070.4 


noise 


AFOE95 


2.5 
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Table 3: Oscillation frequencies for r\ Boo (/xHz) 







I = 





I = 1 


I = 2 


n = 


13 








608.1 (0.4) 


n = 


14 


611.0 


(0.5) 


629.4 (0.3) 




n = 


15 


651.2 


(0.6) 


670.1 (0.5) 




n = 


16 


690.8 


(0.6) 


711.8 (0.4) 


728.4 (0.6) 


n = 


17 


732.6 


(0.4) 


749.3 (0.4)* 












753.4 (0.5)* 




n = 


18 






793.1 (0.7) 


810.5 (0.4) 


n = 


19 


813.1 


(0.4) 




849.9 (0.7) 


n = 


20 


853.6 


(0.3) 






n = 


21 


894.2 


(0.6) 






n = 


22 






955.6 (0.8) 


971.7 (0.6) 


n = 


23 


974.5 


(0.7) 






n = 


24 






1034.3 (0.7) 




A^/ 


40.45 (0.07) 


40.89 (0.19) 


40.41 (0.10) 



Table 4: Frequency separations for rj Boo 



Az; 40.47 ± 0.05 /zHz 

(5^02 3.00 ± 0.35 (Mz 

Sv i 0.78 ± 0.45 (Mz 

D 0.49 ± 0.06 (Mz 

e 1.09 ±0.02 



